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Abstract 

The  problem  of  approximating  convex  bodies  by  polytopes  is  an  important  and  well  studied  prob¬ 
lem.  Given  a  convex  body  K  in  Rd ,  the  objective  is  to  minimize  the  number  of  vertices  (alternatively, 
the  number  of  facets)  of  an  approximating  polytope  for  a  given  Hausdorff  error  e.  Results  to  date 
have  been  of  two  types.  The  first  type  assumes  that  K  is  smooth,  and  bounds  hold  in  the  limit  as 
e  tends  to  zero.  The  second  type  requires  no  such  assumptions.  The  latter  type  includes  the  well 
known  results  of  Dudley  (1974)  and  Bronshteyn  and  Ivanov  (1976),  which  show  that  in  spaces  of 
fixed  dimension,  O^dianJAJ/e)^-1-1/2)  vertices  (alt.,  facets)  suffice.  Our  results  are  of  this  latter 
type. 

In  our  first  result,  under  the  assumption  that  the  width  of  the  body  in  any  direction  is  at  least 
e,  we  strengthen  the  above  bound  to  0(Uarea(A')/e^d_1^2).  This  is  never  worse  than  the  previous 
bound  (by  more  than  logarithmic  factors)  and  may  be  significantly  better  for  skinny  bodies.  Our 
analysis  exploits  an  interesting  analogy  with  a  classical  concept  from  the  theory  of  convexity,  called 
the  Mahler  volume.  This  is  a  dimensionless  quantity  that  involves  the  product  of  the  volumes  of  a 
convex  body  and  its  polar  dual. 

In  our  second  result,  we  apply  the  same  machinery  to  improve  upon  the  best  known  bounds  for 
answering  ^-approximate  polytope  membership  queries.  Given  a  convex  polytope  P  defined  as  the 
intersection  of  halfspaces,  such  a  query  determines  whether  a  query  point  q  lies  inside  or  outside  P, 
but  may  return  either  answer  if  q' s  distance  from  P’s  boundary  is  at  most  e.  We  show  that,  without 
increasing  storage,  it  is  possible  to  reduce  the  best  known  search  times  for  e-approximate  polytope 
membership  significantly.  This  further  implies  improvements  to  the  best  known  search  times  for 
approximate  nearest  neighbor  searching  in  spaces  of  fixed  dimension. 


1  Introduction 

Approximating  convex  bodies  by  polytopes  is  a  fundamental  problem,  which  has  been  extensively  studied 
in  the  literature.  (See  (6)  for  a  recent  survey.)  The  problem  is  to  determine  the  minimum  number  of 
vertices  n  (alternatively,  the  minimum  number  of  facets)  of  an  approximating  polytope  for  a  given  error 
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e  >  0.  Error  is  commonly  measured  through  the  Hausdorff  distance  (6).  Given  a  convex  body  K  in 
d-dimensional  space,  we  say  that  a  polytope  P  e- approximates  K  if  the  Hausdorff  distance  between  K 
and  P  is  at  most  e.  Throughout  this  paper,  we  will  restrict  attention  to  the  Hausdorff  metric  and  to 
Euclidean  space  of  constant  dimension  d. 

Although  much  work  has  been  done  on  the  optimization  problem  of  finding  the  “best”  approximating 
polytope  for  a  specific  input  according  to  various  criteria  (see,  for  example,  (9;  13;  17;  22)),  our  interest 
is  in  establishing  bounds  on  the  combinatorial  complexity  of  approximating  general  convex  bodies.  Ap¬ 
proximation  bounds  are  typically  expressed  as  the  number  n  of  vertices  or  facets  of  a  convex  polytope 
needed  to  achieve  a  certain  approximation  error  e  for  a  given  convex  body  K.  Results  are  of  two  common 
types.  In  both  cases,  it  is  shown  that  there  exists  £o  >  0,  and  the  bounds  hold  for  all  £  <  £0.  In  the  first 
type,  which  we  call  nonuniform  bounds,  the  value  of  £o  depends  on  K  (for  example,  on  K' s  maximum 
curvature).  Such  bounds  are  often  stated  as  holding  “in  the  limit”  as  £  approaches  zero,  or  equivalently 
as  n  approaches  infinity.  In  the  second  type,  which  we  call  uniform  bounds,  the  value  of  £o  is  independent 
of  K. 

As  an  example  of  a  nonuniform  bound,  consider  the  following  result  due  to  Gruber  (14).  Let  dK 
denote  K’s  boundary.  For  a  point  x  €  dK,  let  k(x)  denote  the  Gaussian  curvature  at  x,  and  let  da{x) 
denote  the  ordinary  differential  surface  element  at  x.  Under  certain  smoothness  conditions  on  K ,  Gruber 
showed  that  there  exists  an  approximating  polytope  with  n  vertices,  where 


as  £  — ■>  0.  The  constant  hidden  in  the  O-notation  depends  only  on  dimension.  This  bound  also  applies 
to  the  number  of  facets  n  in  an  approximating  polytope.  Recently,  Clarkson  (10)  has  extended  these 
results  and  shown  a  similar  bound  on  the  triangulation  complexity  for  more  general  manifolds  that  may 
not  represent  the  boundary  of  a  convex  body.  Other  examples  of  nonuniform  bounds  can  be  found  in 
(15;  21;  26;  27). 

In  contrast,  examples  of  uniform  bounds  include  Dudley  (12)  and  Bronshteyn  and  Ivanov  (5).  These 
bounds  hold  without  any  smoothness  assumptions,  and  they  hold  even  for  polytopes.  Dudley  showed 
that,  for  £  <  1,  any  convex  body  K  can  be  £-approximated  by  a  polytope  P  with  0((diam(A')/£)(d_1)/2) 
facets.  Bronshteyn  and  Ivanov  showed  the  same  bound  holds  for  the  number  of  vertices.  These  bounds 
are  tight  up  to  constant  factors  (e.g.,  when  K  is  a  ball).  Again,  constants  hidden  in  the  O-notation 
depend  only  on  d.  These  results  have  many  applications,  for  example,  in  the  construction  of  coresets  (1). 
In  this  paper,  we  consider  only  uniform  bounds. 

1.1  Our  Results 

We  present  two  main  results,  which  arise  by  replacing  the  approximation  bounds  of  both  Dudley  and 
Bronshteyn  and  Ivanov  with  bounds  that  are  more  sensitive  to  the  local  geometry  of  the  convex  body. 
Our  first  result  shows  that  the  existing  bounds  can  be  significantly  improved  for  the  case  when  I\  has 
directional  width  smaller  than  the  diameter  along  some  directions,  that  is,  when  K  is  skinny.  It  seems 
intuitive  that  skinny  bodies  should  be  easier  to  approximate,  but  we  are  unaware  of  any  uniform  bounds 
in  the  literature  that  capture  this  fact.  Under  the  reasonable  assumption1  that  the  width  of  the  body 
in  any  direction  is  at  least  £,  we  show  that  it  is  possible  to  construct  an  ^approximating  polytope  P 
in  which  the  number  of  facets2  is  0(y/ area(A') / £(rf~1)/2) ,  where  area(AT)  denotes  the  surface  area  of  K. 
Although  we  do  not  include  it  here,  a  similar  bound  can  be  proved  for  the  number  of  vertices.  Since 
area(.A)  is  bounded  by  0(diam(AT)d_1),  our  bound  is  as  good  as  Dudley’s  for  fat  convex  bodies  (up  to 
logarithmic  factors),  and  it  is  superior  for  skinny  ones.  These  results  are  presented  in  Section  4,  with 
the  main  result  given  in  Theorem  4.1. 

Our  second  result  involves  approximate  polytope  membership  queries.  Given  a  convex  polytope 
P,  defined  as  the  intersection  of  halfspaces,  an  £-approximate  polytope  membership  query  determines 

1  Without  this  assumption  it  is  impossible  to  achieve  area-based  bounds.  For  example,  given  a  fat  d-dimensional  convex 
body,  embed  it  in  dimension  d  +  2  and  take  its  Minkowski  sum  with  a  ball  of  radius  8  -C  e.  The  surface  area  of  the  resulting 
body  can  be  made  arbitrarily  small  as  a  function  of  5,  but  clearly  the  complexity  of  any  approximating  body  is  not  smaller 
than  for  the  original  body.  Conversely,  if  the  body  is  of  width  less  than  e  in  some  direction,  then  by  projecting  the  body 
onto  a  hyperplane  orthogonal  to  this  direction,  we  can  reduce  the  problem  to  a  convex  approximation  problem  in  one  lower 
dimension.  This  can  be  repeated  until  the  body’s  width  is  sufficiently  large,  after  which  our  bounds  apply. 

2Throughout,  we  treat  e  as  an  asymptotic  quantity  and  use  O  to  indicate  that  logarithmic  factors  in  e  are  ignored. 
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whether  a  query  point  q  lies  inside  or  outside  P,  but  it  may  return  either  answer  if  q's  distance  from  P’s 
boundary  is  at  most  e.  Polytope  membership  queries  find  applications  in  many  geometric  areas,  such 
as  linear  programming  queries,  ray  shooting,  nearest  neighbor  searching,  and  the  computation  of  convex 
hulls  (7;  8;  18;  20;  23).  In  dimension  d  <  3,  it  is  possible  to  build  a  data  structure  of  linear  size  that  can 
answer  such  queries  in  logarithmic  time  (11).  In  higher  dimensions,  however,  all  exact  data  structures 
with  roughly  linear  space  take  0(n1_1/L<V2J )  query  time  (19),  which,  except  in  small  dimensions,  is  little 
better  than  brute-force  search. 

Assuming  that  P  has  unit  diameter,  a  straightforward  application  of  Dudley’s  method  allows  such 
queries  to  be  answered  in  0(l/el'd~1^2)  storage  and  query  time  (by  a  brute-force  enumeration  of  the 
approximating  halfspaces).  We  recently  showed  that  there  is  a  very  simple  hierarchical  search  algorithm 
that  achieves  the  same  storage  bound  but  improves  the  query  time  to  0(  1/e^-1^4),  thus  reducing  the 
exponent  by  half  (2).  In  this  paper,  we  show  that,  with  the  same  storage  bounds,  it  is  possible  to  again 
reduce  the  exponent  in  the  query  time  by  half,  to  0(l/e(d_1^8).  We  also  present  a  space-time  tradeoff 
that  improves  query  times  to  the  square  root  of  the  previous  bounds  throughout  the  spectrum.  The 
results  are  summarized  in  Figure  1(a). 


Figure  1:  The  multiplicative  factor  in  the  exponent  of  the  1/e  term  for  (a)  polytope  membership  queries 
and  (b)  approximate  nearest  neighbor  (ANN)  queries.  The  0(1)  term  in  the  exponent  corresponds  to  a 
constant  that  does  not  depend  on  d. 

It  is  notable  that  the  improvement  involves  no  change  to  the  simple  data  structure  of  (2),  just  a 
better  analysis  of  its  complexity.  These  results  are  presented  in  Section  5,  and  the  main  results  is  given 
in  Theorem  5.1.  As  shown  in  (2),  these  results  provide  immediate  improvements  to  the  best  space-time 
tradeoffs  for  approximate  nearest  neighbor  searching  (see  Figure  1(b)). 

1.2  Overview  of  Methods 

Our  results  stem  from  a  locally  sensitive  approach  to  convex  approximation.  Our  area-based  bound  on 
the  complexity  of  convex  body  approximation  comes  about  through  a  consideration  of  the  complexity 
of  approximating  small  patches  on  the  boundary  of  the  convex  body,  called  caps.  Given  a  halfspace 
H  whose  bounding  hyperplane  intersects  K ,  we  define  a  cap  to  be  the  portion  of  dK  that  is  contained 
within  H .  The  width  of  the  cap  is  the  maximum  distance  between  any  point  on  the  cap  and  the  bounding 
hyperplane  of  H ,  and  an  e-cap  is  a  cap  of  width  e.  A  set  of  points  is  said  to  stab  all  £-caps  if  every 
e-cap  has  a  nonempty  intersection  with  the  set.  Bronshteyn  and  Ivanov  (5)  observed  that  the  convex 
hull  of  any  set  of  points  that  stabs  all  £-caps  yields  an  ^approximation  of  K.  Assume  that  K  has  unit 
diameter,  and  let  S  be  the  (d  —  l)-dimensional  sphere  of  radius  3  centered  at  the  origin,  which  we  call 
the  Dudley  hypersphere.  (Any  sufficiently  large  constant  radius  suffices.)  We  may  associate  any  £-cap  C 
with  the  surface  patch  on  S  consisting  of  points  whose  closest  point  on  dK  lies  on  C;  we  call  this  the 
Voronoi  patch  of  C. 
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Our  key  geometric  observation  is  that  the  product  of  the  area  of  any  £-cap  and  its  associated  Voronoi 
patch  must  be  large.  Intuitively,  if  the  surface  area  of  the  £-cap  is  small,  then  its  curvature  must  be 
high,  and  so  its  associated  Voronoi  patch  must  have  relatively  large  area.  We  exploit  this  observation  to 
efficiently  stab  all  the  £-caps  using  a  two-pronged  strategy.  We  stab  £-caps  that  have  large  surface  area 
by  selecting  points  randomly  on  the  boundary  of  K  using  density  proportional  to  the  surface  area.  We 
stab  the  £-caps  that  have  small  surface  area  by  selecting  points  randomly  on  the  Dudley  hypersphere 
(again  using  density  proportional  to  the  surface  area),  and  then  choosing  their  nearest  neighbor  on  dK. 
We  will  see  in  Section  4  that  this  allows  us  to  stab  all  the  £-caps  with  0(i/area(A')/£^”1)/2)  points. 

Our  approach  to  answering  polytope  membership  queries  for  a  convex  body  K  is  based  on  a  quadtree 
search,  where  the  basic  problem  involves  approximating  the  intersection  of  the  boundary  of  K  with  a 
single  quadtree  box.  Our  improved  bounds  for  polytope  membership  queries  follow  as  a  consequence  of 
establishing  improved  bounds  for  such  approximations,  whenever  the  diameter  of  the  quadtree  box  is  at 
least  yfe.  Roughly  speaking,  we  show  that  if  A\  denotes  the  area  of  dK  (~l  b,  and  A 2  denotes  the  area 
of  the  associated  Voronoi  patch,  then  dK  fl  b  can  be  £-approximatecl  using  0(y/ A\Ail e^-1)/2)  facets. 
From  this  it  is  not  hard  to  show  that  for  a  box  b  of  diameter  yfe,  if  Dudley’s  method  uses  m  facets  to 
£-approximate  dK  (~l  b,  then  0(y/m)  facets  suffice.  We  then  exploit  this  fact  in  presenting  an  improved 
analysis  of  the  search  algorithm  of  (2). 

The  geometric  underpinnings  needed  to  show  that  the  product  of  areas  is  bounded  is  presented  in 
Section  3.  We  introduce  a  dual  notion  of  cap,  called  a  dual  cap,  and  we  explore  the  relationship  between 
the  areas  of  surface  caps  and  their  Voronoi  patches.  We  demonstrate  an  interesting  parallel  between 
these  concepts  and  a  classical  notion  from  the  theory  of  convexity,  called  the  Mahler  volume  (4;  25). 
The  Mahler  volume  of  a  convex  body  is  a  dimensionless  quantity  that  involves  the  product  of  the  body’s 
volume  and  the  volume  of  its  polar  body.  We  demonstrate  that  caps  and  their  Voronoi  patches  exhibit 
a  similar  polar  relationship. 

The  improvements  presented  here  do  not  arise  simply  by  a  more  careful  application  of  the  analysis 
given  in  (2).  The  analysis  given  there  relied  on  a  direct  application  of  Dudley’s  and  Bronshteyn  and 
Ivanov’s  bounds.  The  key  to  the  improvements  of  this  paper  is  the  development  of  uniform  bounds  that 
are  sensitive  to  local  geometry.  Sensitivity  to  local  geometry  has  been  applied  in  asymptotic  contexts 
(under  the  assumption  of  smoothness)  by  Gruber  (14)  and  Clarkson  (10).  It  has  also  been  applied  for 
uniform  bounds,  but  only  in  the  two-dimensional  setting  by  Rote  (24).  The  techniques  developed  here 
provide  powerful  tools  for  analyzing  the  complexity  of  approximating  arbitrary  convex  bodies  on  a  local 
level,  and  we  believe  that  these  techniques  will  be  of  wider  applicability. 

2  Preliminaries 

Let  K  denote  a  convex  body  in  that  is,  a  nonempty,  closed,  bounded  convex  set.  Let  dK  denote 
K's  boundary.  Given  6  >  0,  we  say  that  a  convex  body  is  5-smooth  if  at  each  point  q  £  dK  there  exists 
a  ball  of  radius  6  that  lies  entirely  within  K  and  contains  q  on  its  boundary.  We  say  that  K  is  smooth 
if  it  is  (5-smooth  for  some  6  >  0.  If  K  is  smooth,  then  for  any  q  £  dK,  there  is  a  unique  supporting 
hyperplane  for  K  containing  q,  which  we  denote  by  h(q). 

Given  any  bounded  set  X  in  Rd,  which  for  us  will  be  a  convex  body,  we  use  vol(A)  to  denote  its 
d-dimensional  Hausdorff  measure.  Given  a  (d—  l)-dimensional  manifold  Y  in  Rd,  for  example,  a  surface 
patch  on  a  convex  body,  we  let  area(V)  denote  its  (d—  l)-dimensional  Hausdorff  measure.  Given  a  convex 
body  X  in  Rd,  we  use  area(X)  as  a  shorthand  for  area(cGT).  We  say  that  a  polytope  P  e- approximates 
K  if  the  Hausdorff  distance  between  K  and  P  is  at  most  s. 

Let  H  be  a  halfspace  whose  bounding  hyperplane  intersects  K.  The  intersection  K  n  H  is  called  the 
cap  generated  by  H .  Since  we  are  interested  in  the  surface  of  K,  henceforth  we  will  use  the  term  cap 
in  a  nonstandard  sense  as  the  portion  of  the  boundary  of  the  cap  that  coincides  with  dK  (in  essence, 
the  outer  skin  of  the  cap).  Let  q  denote  a  point  of  the  cap  of  maximum  distance  from  dH.  By  local 
maximality,  h{q)  and  dH  are  parallel  to  each  other,  and  the  distance  between  them  is  called  the  width 
of  the  cap.  A  cap  of  width  £  is  called  an  £-cap.  We  refer  to  q  as  the  cap’s  apex,  and  K  n  dH  as  the  cap’s 
base,  which  we  denote  by  T?.  By  smoothness,  for  each  q  £  dK,  there  is  a  unique  £-cap  having  q  as  its 
apex,  which  we  denote  by  Cq  (see  Fig  2(a)). 

Given  any  point  q  external  to  K,  the  dual  cap3  generated  by  q  is  the  topological  closure  of  the 

3The  term  “dual  cap”  arises  from  the  observation  that,  given  an  inclusion-reversing  dual  transformation,  if  we  let  K* 
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Figure  2:  Caps  and  dual  caps. 


intersection  of  all  halfspaces  that  contain  K  but  not  q.  As  before,  since  we  are  interested  in  the  surface 
of  K,  henceforth  we  use  the  term  dual  cap ,  denoted  by  Dq,  to  refer  to  the  portion  of  the  boundary  of  the 
dual  cap  that  coincides  with  dK  (which  can  be  visualized  as  the  subset  of  points  of  dK  illuminated  by 
a  light  source  at  q).  Let  q'  be  the  closest  point  of  K  to  q  (see  Fig  2(b)). 

The  distance  from  q  to  q'  is  called  the  width  of  the  dual  cap.  A  dual  cap  of  width  £  is  called  an 
£-dual  cap.  The  point  q  is  the  dual  cap’s  apex,  and  h(q')  fl  conv({g}  U  K)  is  the  dual  cap’s  base,  which 
we  denote  by  Tq. 


Figure  3:  Voronoi  cells:  (a)  associated  with  a  surface  patch  X  and  (b)  a  dual  cap  D  and  the  associated 
Voronoi  cell  on  H(q). 

Given  a  surface  patch  X  C  dK,  let  Vor(X)  denote  the  set  of  points  exterior  to  K  whose  closest  point 
on  dK  lies  within  X.  Let  S  denote  the  (d—  l)-dimensional  sphere  of  radius  3  centered  at  the  origin,  which 
we  call  the  Dudley  hyper  sphere.  (The  name  is  inspired  by  the  sphere  used  in  polytope  approximation 
construction  (12).  The  value  3  can  be  replaced  by  any  constant  that  is  at  least  a  constant  factor  larger 
than  1  +  £.)  Throughout,  we  assume  that  K  has  been  scaled  and  translated  to  lie  within  the  ball  of  unit 
radius  centered  at  the  origin.  It  follows  that  the  distance  from  any  point  of  S  to  its  closest  point  of  K  is 
at  least  2.  Given  a  cap  or  dual  cap  X,  Vor(A)  n  5  is  the  surface  patch  consisting  of  points  of  S  whose 
closest  point  on  dK  lies  within  X  (see  Figure  3(a)). 

Because  we  are  interested  in  the  local  geometry  of  K ’s  boundary,  it  will  be  useful  to  limit  the  diameter 
of  caps  and  dual  caps.  Given  S  >  0,  define  Bg(q)  to  be  a  closed  Euclidean  ball  of  radius  S  centered  at 
q.  Define  the  5 -restriction  of  Cq,  denoted  Cq  ,  to  be  Cq  D  Bg(q).  Similarly,  define  the  8-restriction  of  Dq, 
denoted  Dq  ,  to  be  Dq  fl  Bg(q). 

3  Caps,  Dual  Caps,  and  Areas 

In  this  section  we  present  the  principal  geometric  underpinnings  of  our  approximation  results.  In  par¬ 
ticular,  we  consider  the  relationship  between  the  surface  areas  of  both  restricted  and  unrestricted  £-dual 
caps  and  the  surface  areas  of  their  associated  Voronoi  patches  on  S.  As  mentioned  in  the  introduction, 
an  £-dual  cap  of  small  area  must  have  high  curvature,  and  hence  its  Voronoi  patch  on  the  Dudley  hyper¬ 
sphere  should  be  large.  We  will  establish  this  relationship  for  dual  caps  by  providing  a  lower  bound  on 
the  product  of  the  areas  of  the  £-dual  cap  Dq  and  its  associated  Voronoi  patch,  Vor (D  )  fl  S.  Analogous 

and  q*  be  the  respective  duals  of  K  and  q ,  the  dual  cap  of  K  generated  by  q  is  the  dual  of  the  cap  of  K*  generated  by  q* . 
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results  can  be  shown  for  caps,  but  dual  caps  suffice  for  our  purposes.  Our  analysis  will  reveal  an  interest¬ 
ing  similarity  between  these  products  and  the  Mahler  volume  of  a  convex  body  (4;  16;  25).  In  Section  3.1 
we  start  with  the  simpler  unrestricted  case,  and  in  Section  3.2  we  extend  this  to  the  restricted  case. 

3.1  Area  Product  Bounds:  Unrestricted  Case 

In  this  section  we  limit  attention  to  unrestricted  dual  caps.  Rather  than  dealing  with  spherical  patches, 
it  will  be  easier  to  start  with  hyperplane  patches  and  then  generalize  to  the  spherical  case.  Given  a  point 
q  external  to  K,  let  q'  be  the  closest  point  on  dK  to  q ,  and  recall  that  h(q')  is  the  unique  supporting 
hyperplane  at  q' .  Let  H(q)  be  a  hyperplane  that  is  parallel  to  h(q'),  tangent  to  the  Dudley  hypersphere 
S,  and  on  the  same  side  of  K  as  h(q')  (see  Figure  3(b)).  We  call  this  the  Dudley  hyperplane  for  q.  We 
prove  two  bounds  regarding  the  product  of  the  area  of  dual  caps  and  their  Voronoi  patches,  one  for  the 
Dudley  hyperplane  and  one  for  the  Dudley  hypersphere. 

Lemma  3.1  Let  0  <  e  <  1.  Let  I\  be  a  smooth,  convex  body  that  lies  within  a  unit  ball  in  M  ,  and  let  q 
be  a  point  at  distance  e  from  I\.  Let  Dq  denote  the  e-dual  cap  generated  by  q.  Then  there  is  a  constant 
c,  depending  only  on  dimension,  such  that: 

( i )  area (D  )  ■  area(Vor (D  )  H  H(q))  >  c  ■  ed~x ,  where  H(q)  is  the  Dudley  hyperplane  for  q, 

(ii)  if  K  is  e-smooth,  then  area (D)  ■  area(Vor(I?<3)  n  S)  >  c  ■  £d_1,  where  S  is  the  Dudley  hypersphere 
for  K . 

Before  considering  Vor (D  )  ft  H(q),  it  will  be  helpful  to  define  a  body  that  is  related  to  the  polar  dual 
of  the  base  Tq.  The  parallel  hyperplanes  h(q')  and  H(q)  can  each  be  identified  with  M“_1,  by  taking  the 
origin  to  be  located  at  the  orthogonal  projection  of  q  onto  each  hyperplane.  (In  the  case  of  h(q')  this  is 
q' .)  The  coordinate  axes  for  both  are  taken  to  be  any  consistent  set  of  d  —  1  orthonormal  vectors  lying 
parallel  to  these  hyperplanes.  Since  T9  is  the  intersection  of  a  (d  —  l)-dimensional  hyperplane  with  a 
bounded  convex  set  that  contains  q'  (the  convex  hull  of  K  U  {?}),  we  may  view  it  as  a  convex  body  in 
Kd_1  that  contains  the  origin.  Let  ||u||  denote  the  Euclidean  length  of  v  £  Kd_1.  For  each  v  £  r?,  we 
define  hv  to  be  a  halfspace  in  containing  the  origin  whose  bounding  hyperplane  is  orthogonal  to  v 
and  at  distance  e/||u||  from  the  origin  (on  the  same  side  of  the  origin  as  v ). 

Define  Tq  =  n„er,  which  we  treat  as  a  subset  of  H(q).  This  body  is  bounded,  convex,  and 

contains  the  origin.  Rather  than  analyzing  Vor (D  )  C\H(q)  directly,  we  will  focus  on  T  .  First,  we  show 
that  r?  is  a  subset  of  Vor(D  )  Cl  H(q). 


(b) 

Figure  4:  Proof  of  Lemma  3.2. 


Lemma  3.2  Given  the  preconditions  of  Lemma  3.1,  r?  C  Vor (Dq)  ft  H(q). 
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Proof :  Consider  any  point  w  £  dK  \Dq.  It  suffices  to  show  that  w  cannot  be  the  closest  point  of  dK  to 
any  point  p  £  Tq.  It  will  simplify  matters  to  consider  the  problem  in  a  2-dimensional  setting.  Consider 
the  unique  plane  $  containing  the  points  q ,  q' ,  and  w.  For  the  sake  of  visualization,  we  identify  this 
plane  with  R2  so  that  q'  is  the  origin  and  the  axes  are  positioned  so  that  q  lies  on  the  y-axis  at  distance  e 
above  q'  (see  Figure  4(a)).  Let  K' ,  h'(q')  and  H'{q)  denote  the  respective  intersections  of  K ,  h(q'),  and 
H{q)  with  this  plane,  and  let  p'  denote  the  orthogonal  projection  of  p  onto  <I>.  Without  loss  of  generality, 
we  may  assume  that  w  lies  in  the  positive- a;  halfplane.  Let  t  be  the  point  of  tangency  of  dK'  with  respect 
to  q  in  this  halfplane.  Let  t'  denote  the  point  at  which  the  line  qt  intersects  h'(q'). 

Clearly,  any  point  of  dK'  can  be  closest  to  some  point  of  r?  only  if  it  lies  on  the  upper  hull  of  K 
(that  is,  the  portion  of  dK'  visible  from  a  point  at  y  =  +00).  Throughout,  we  restrict  attention  to  the 
portion  of  the  upper  hull  of  K'  that  lies  in  the  positive- a:  halfplane. 

Given  a  point  u  £  4>,  let  ux  and  uy  denote  its  coordinates  relative  to  the  above  coordinate  system.  Let 
v  denote  the  horizontal  vector  from  q'  to  t' ,  whose  length  is  clearly  t'x.  By  viewing  v  as  a  vector  in  ]Rd_1 
on  the  plane  h(q'),  it  contributes  a  halfspace  hv  to  the  body  r?,  which  is  bounded  by  a  hyperplane  that 
is  orthogonal  to  v  (and  hence  orthogonal  to  4>)  and  lies  at  distance  e/t'x  from  the  orthogonal  projection 
of  q  onto  H'(q)  (see  Figure  4(b)). 

To  complete  the  proof,  we  will  show  that  w  is  not  the  closest  point  on  dK'  to  p' .  To  see  why  this 
suffices,  observe  that  the  squared  Euclidean  distance  from  p  to  any  point  on  4)  can  be  decomposed  into 
two  components,  one  parallel  to  4>  and  one  orthogonal  to  4>.  The  orthogonal  component  is  the  same  for 
all  points  on  4>,  and  so  may  be  ignored  when  comparing  distances.  The  parallel  component  is  determined 
by  the  distances  to  p' .  Thus,  if  w  is  not  the  closest  point  on  dK'  to  //,  then  w  is  not  the  closest  point 
of  dK  to  p. 

Since  K'  is  convex  and  smooth,  there  is  a  support-based  condition  for  any  point  u  £  dK'  to  be  the 
closest  point  on  dK'  to  a  point  p'  £  K' .  In  particular,  the  line  segment  up’  must  be  perpendicular  to 
the  (unique)  supporting  hyperplane  to  K'  at  u.  To  apply  this,  consider  the  convex  body  conv({q}  U  K'), 
and  let  Tq  denote  the  surface  patch  on  this  boundary  consisting  of  the  union  of  closed  line  segments 
joining  q  to  its  points  of  tangency  on  K' .  (We  just  show  the  right-side  tangent  in  Figure  4(b).)  Based 
on  the  support-based  criterion  given  above,  is  easy  to  see  that,  with  respect  to  the  body  conv({<7}  U  K), 
Vor {Tq)  C\H{q)  is  the  same  as  Vor(D  )  fl  H{q).  (Vor (Tq)  fl  H{q)  is  shown  as  a  solid  line  segment  on  H{q) 
in  Figure  4(b).  Although  Tq  is  not  smooth,  it  can  be  made  so  by  taking  its  Minkowski  sum  with  a  ball 
of  infinitesimally  small  radius.) 

To  show  that  w  is  not  the  closest  point  on  dK’  to  p' ,  we  will  show  that  \\qp'\\  <  ||u;f/||.  To  this 
end,  consider  a  ray  in  4>  emanating  from  q  in  the  direction  of  an  outward  normal  to  the  segment  qt. 
The  absolute  value  of  the  slope  of  the  line  qt  is  e/t'x,  and  therefore,  the  slope  of  this  ray  is  t'x/e.  Let  i 
denote  the  distance  from  q  to  H(q).  Since  each  point  of  K  is  at  distance  at  least  2  from  H(q ),  we  have 
£  >  2  —  £  >  1,  it  follows  that  this  ray  intersects  H'(q)  at  a  distance  of  £s/t’x  from  the  y-axis,  which  is  to 
the  right  of  the  point  at  which  hv’ s  bounding  hyperplane  intersects  H'(q).  Since  p  £  hv,  p'  lies  to  the 
left  of  this  ray.  It  is  easy  to  see,  however,  that,  restricted  to  the  positive-a:  halfplane,  any  point  lying  to 
the  left  of  the  ray  is  closer  to  q  than  to  any  point  lying  below  the  line  qt.  Since  w  £  Dq ,  it  lies  below  qt, 
and  therefore,  we  achieve  the  desired  conclusion  that  w  is  not  the  closest  point  on  dK’  to  p' .  □ 

In  order  to  bound  the  product  of  the  areas  of  Tq  and  r?,  we  will  demonstrate  the  relationship  between 
r?  and  the  polar  dual  of  T  ,  and  then  apply  a  well  known  result  from  the  theory  of  convexity,  called 
the  Mahler  volume.  First,  we  recall  some  facts  about  the  polar  transform.  Given  a  vector  v  £  Rd  other 
than  the  origin,  define  polar(u)  to  be  the  halfspace  that  contains  the  origin  whose  bounding  hyperplane 
is  orthogonal  to  v  and  at  distance  l/||u||  from  the  origin  (on  the  same  side  of  the  origin  as  v).  Generally, 
given  r  >  0,  define  polarr(w)  to  be  the  halfspace  at  distance  r2/||i;||  (see  Figure  5(a)).  Given  a  convex 
body  K  that  contains  the  origin  in  its  interior,  define  polarr(iv)  to  be  the  intersection  of  the  halfspaces 
polarr(u),  for  all  v  £  I\  (see  Figure  5(b)).  Clearly,  polarr(AT)  is  a  scaled  copy  of  polar(A')  by  a  factor  of 
r2. 

The  bodies  polar  j£{K)  and  r?  are  bounded  by  the  same  set  of  halfspaces,  and  so  we  have: 

Lemma  3.3  Given  a  smooth  convex  body  K  and  an  exterior  point  q,  T?  =  polar ^(rg). 

Given  the  above  results,  we  can  now  prove  the  first  part  of  Lemma  3.1.  Recall  that  0  <  e  <  1. 
Since  r?  is  a  subset  of  the  orthogonal  projection  of  Dq  onto  the  plane  h(q'),  area(T9)  <  area(Z?9).  By 
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Figure  5:  The  generalized  polar  transform  and  polar  body. 


Lemmas  3.2  and  3.3,  we  have 

area(Vor(.D  )  fl  H(q))  >  area(r?)  =  area(polarv^(rg)). 

As  observed  earlier,  polarr(r?)  is  a  scaled  copy  of  polar(T9)  by  a  factor  of  r2.  Since  this  is  a  (d  —  1)- 
dimensional  set,  we  have  area(polar ^(Tg))  =  £d_1 -area(polar(r9)).  By  known  inequalities  on  the  Mahler 
volume  (16),  there  is  a  constant  c  depending  only  on  d  such  that  area(T9)-area(polar(T9))  >  c.  Combining 
these  observations,  we  have 

area(-Dg)  •  area(Vor(.D(J)  fl  H(q)) 

>  area(r?)  •  (ed_1  •  area(polar(T?))) 

>  C  •  £d_1, 

which  establishes  Lemma  3.1(i). 

In  order  to  establish  the  second  part  of  Lemma  3.1,  we  need  to  assume  that  K  is  sufficiently  smooth. 
If  not,  r  might  have  very  small  area  (much  smaller  than  e d~1),  but  the  entire  surface  area  of  the  Dudley 
hypersphere  is  bounded  above  by  a  constant.  If  K  is  £-smooth,  then  it  is  easy  to  show  that  Tg  contains 
a  (d  —  l)-dimensional  ball  (centered  at  q')  of  radius  at  least  e/\/Z. 

Consider  the  area  of  T  .  Each  of  its  bounding  halfspaces  hv  is  defined  by  a  hyperplane  whose  distance 
from  the  origin  is  £/||u||.  By  smoothness,  along  any  direction,  there  exists  a  vector  v  £  T  of  length  at 
least  e/y/3,  so  the  associated  bounding  hyperplane  hv  for  T?  is  at  distance  at  most  £/(e/V3)  =  a/3  from 
the  origin  on  H(q).  The  projection  of  any  differential  element  on  H(q)  lying  within  a  constant  distance 
of  the  origin  undergoes  a  constant  factor  decrease4  when  projected  onto  to  the  Dudley  hypersphere. 
Thus,  for  some  constant  factor  d ,  we  have  area(Vor(DlJ)  D  S)  >  d  ■  area(Vor(.D  )  fl  H{q)).  Combining 
this  with  Lemma  3.1(i)  establishes  Lemma  3.1  (ii) . 

3.2  Area  Product  Bounds:  Restricted  Case 

In  this  section  we  generalize  the  results  of  the  previous  section  to  the  case  of  restricted  dual  caps. 
(As  before,  analogous  results  hold  for  restricted  caps,  but  these  results  are  not  needed  here.)  Recall 
that,  given  a  parameter  6  >  0,  the  ^-restricted  cap  Cd  and  the  (5-restrictecl  dual  cap  come  about  by 
intersecting  the  standard  cap  and  dual  cap,  respectively,  with  a  ball  of  radius  <5  centered  at  q.  Restriction 
complicates  the  simple  relationship  between  T?  and  the  polar  body.  Although  the  analysis  of  the  section 
is  structurally  similar  to  that  of  the  previous  section,  there  are  added  technical  complications. 

As  before,  we  start  by  considering  the  intersection  of  the  Voronoi  region  with  the  Dudley  hyperplane 
H(q)  and  then  generalize  to  the  Dudley  hypersphere  S. 

Lemma  3.4  Let  0  <  £  <  1  and  3 e  <  6  <  3 yfe.  Let  I\  be  a  smooth,  convex  body  that  lies  within  a  unit 
ball  in  and  let  q  be  a  point  at  distance  e  from  I\.  Let  D 5  denote  the  5-restriction  of  the  e-dual  cap 
generated  by  q.  Then  there  is  a  constant  c,  depending  only  on  dimension,  such  that: 

4Consider  a  differential  element  dv  on  H(q).  The  corresponding  vector  with  respect  to  the  center  of  the  Dudley 
hypersphere  is  of  length  x  =  y/9  +  ||n||2,  and  it  is  at  angle  9  =  arctan(||n||/3)  with  respect  to  vertical.  Projecting  this 
element  onto  the  Dudley  hypersphere  results  in  a  scaling  by  (3/a;)  cos 9.  Since  ||u||  <  \/3,  we  have  x  <  2v3  and  cos 9  >  v3/2. 
Therefore,  the  projection  of  this  differential  element  onto  the  Dudley  sphere  is  at  least  (3/4 )dv,  which  implies  that  the 
projection  results  in  a  constant-factor  decrease  in  the  surface  area. 
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Figure  6:  The  definitions  of  and  F^.  (The  figure  is  not  drawn  to  scale.) 


(z)  area(Dq)  •  area(Vor(.D5)  n  H(q))  >  c  ■  5 2<'d  1\  where  H(q)  is  the  Dudley  hyperplane  for  q, 

(ii)  if  K  is  e-smooth,  then  area (Dq)  ■  area(/Vor(Dq)  Pi  S)  >  c-52^d_1^ ,  where  S  is  the  Dudley  hypersphere 
for  K . 

Define  Tq  to  be  the  ^-restricted  base  of  the  dual  cap,  by  which  we  mean  the  intersection  of  the  base 
with  a  ball  of  radius  S  centered  at  q.  Because  q  lies  at  distance  e  above  h'(q),  T*  is  the  intersection  of 
Tq  and  a  (d  —  l)-dimensional  ball  of  radius  S'  =  V52  —  e2,  centered  at  q'  (see  Figure  6(a)).  As  before, 
we  define  an  intermediate  convex  body  upon  which  to  base  our  analysis.  Recall  that  h{q')  and  H(q) 
can  each  be  identified  with  Kd_1  by  taking  the  origin  to  be  located  at  the  orthogonal  projection  of  q 
onto  each  hyperplane.  We  may  therefore  view  as  a  convex  body  in  ]Rd_1  that  contains  the  origin. 
For  each  v  £  rq,  we  define  a  halfspace  hv  in  Kd_1  as  follows.  If  ||u||  >  S/3,  then  hv  is  the  halfspace 
containing  the  origin  whose  bounding  hyperplane  is  orthogonal  to  v  and  at  distance  5/6  from  the  origin 
(see  Figure  6(b)).  If  ||u||  <  5/3,  then  hv  is  the  (d—  l)-dimensional  halfspace  containing  the  origin  whose 
bounding  hyperplane  is  orthogonal  to  v  and  at  distance  e/(2||u||)  from  the  origin  (see  Figure  6(c)). 

Define  T/  =  n„er«  h",  which  we  treat  as  a  subset  of  H(q).  Note  that,  when  ||i>|  <  5/3,  the 
hyperplane  bounding  hv  in  this  definition  is  closer  to  the  origin  than  that  given  in  the  definition  of 
Tq  from  the  previous  section.  This  body  is  bounded,  convex,  and  contains  the  origin.  Analogous  to 
Lemma  3.2  we  have  the  following. 

Lemma  3.5  Given  the  preconditions  of  Lemma  3-4,  Lq  C  Vor(Dq)  fl  H(q). 

Proof:  Consider  any  point  w  £  dK  \  Dq  .  If  w  ft  D  ,  then  the  proof  follows  from  the  same  argument 
given  in  Lemma  3.2  (with  a  suitable  adjustment  of  constants).  Thus,  it  suffices  to  consider  only  the  case 
where  w  £  D  \  Bs(q). 

As  in  the  earlier  proof,  define  the  plane  $  containing  the  points  q,  q' ,  and  w ,  and  define  K' ,  h'(q'), 
H'(q),  t,  and  t!  as  done  there.  As  before,  it  suffices  to  show  that  w  is  not  the  closest  point  on  K'  to  the 
orthogonal  projection  of  p  onto  <I>,  which  we  denote  by  p' .  Let  v  denote  the  horizontal  vector  from  q'  to 
t! .  By  viewing  v  as  a  vector  in  Kd_1  on  the  plane  h(q'),  it  contributes  a  halfspace  hv  to  the  body  F(y, 
which  is  bounded  by  a  hyperplane  that  is  orthogonal  to  v  (and  hence  orthogonal  to  $). 

We  consider  two  cases,  depending  on  t'x,  the  length  of  v.  In  the  first  case,  assume  that  t'x  >  5/3  (see 
Figure  7(a)).  We  assert  that  wx  >  5/3.  As  mentioned  earlier,  we  may  assume  that  w  ^  Bs(q).  Let  crw 
and  at  denote  the  absolute  values  of  the  slopes  of  the  lines  q'w  and  qt,  respectively.  Since  w  £  Dq  and 
K'  is  convex,  the  line  q'w  must  intersect  the  line  qt,  which  implies  that  aw  <  tr*.  Since  at  =  e/t!x,  and 
t'x  >  5/3  >  e,  we  have  aw  <  1.  Since  w  £  B$(q),  we  have  wx  >  5/V 2  >  <5/3,  as  desired. 

Having  established  that  wx  >  5/3,  consider  the  vector  v.  Since  it  is  of  length  t'x  >  5/3,  Tf  is  a  subset 
of  the  halfspace  hv  (in  H(q))  whose  bounding  hyperplane  is  orthogonal  to  v  and  at  distance  5/6  from  the 
origin  (the  orthogonal  projection  of  q  onto  H(q)).  Since  p'  £  rf,  it  follows  that  it  lies  within  distance 
5/6  of  the  y-axis.  In  contrast,  w  lies  at  more  than  twice  this  distance  from  the  y-axis.  It  follows  directly 
that  p'  is  closer  to  q'  than  to  w,  and  thus,  w  cannot  be  the  closest  point  to  p'  on  dK' . 

In  the  second  case,  assume  that  t'x  <  5/3.  Let  t"  be  the  intersection  of  the  line  qt  with  dBs(q),  and  let 
q"  be  the  point  on  the  y-axis  horizontal  from  t"  (see  Figure  7(b)).  Let  a"  be  the  absolute  slope  of  the  line 
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(£  +  e)e/2t'x 


Figure  7:  Proof  of  Lemma  3.5. 


q't".  By  convexity,  and  the  fact  that  w  (f  Bg(q),  it  follows  that  w  lies  below  the  line  q’t" .  We  assert  that 
a"  >  ut/ 2.  It  is  easy  to  see  that  this  is  equivalent  to  the  inequality  q'y  —  q "  >  (qy  —  q”)/ 2,  or  equivalently 
Qy  ~  Qy  <  ( qy  —  q'y)/ 2.  By  similar  triangles,  this  latter  inequality  is  equivalent  to  \\t'  —  <7||  <  \\t"  —  q\\/2. 
Clearly,  \\t"  —  q||  =  5.  Since  3e  <  5  and  t’x  <  5/ 3,  we  have 


s  r-«  ii 


2  2 

This  establishes  that  a"  >  at/ 2,  as  desired. 

Observe  that  at  =  e/t'x.  Given  that  a"  >  ert/2  —  s/(2tx),  consider  a  ray  emanating  from  q '  that  is 
orthogonal  to  q't".  Letting  £  denote  the  distance  from  q  to  H(q),  by  similar  triangles  we  see  that  this 
ray  intersects  H'{q)  at  a  distance  from  the  y-axis  of  {£  +  e)a"  >  (£  +  e)e/(2t'x).  As  observed  in  the  proof 
of  Lemma  3.2,  £  >  1,  and  so  this  ray  intersects  H'(q)  to  the  right  of  the  point  where  the  hyperplane 
defining  hv  intersects  H'(q ),  whose  distance  from  the  y-axis  is  £/(2||u||)  =  e/{2 t'x).  Since  p  G  it  lies 
to  the  left  of  this  point.  It  is  easy  to  see,  however,  that,  restricted  to  the  positive-a;  halfplane,  any  point 
lying  to  the  left  of  the  ray  is  closer  to  q'  than  to  any  point  that  lies  below  the  line  q't".  Since  w  is  such 
a  point,  we  achieve  the  desired  conclusion  that  w  is  not  the  closest  point  on  dK'  to  p' .  □ 

In  order  to  bound  the  product  of  the  areas  of  and  we  proceed  as  in  the  preceding  section  by 
demonstrating  the  relationship  between  F^  and  the  polar  dual  of  . 

Lemma  3.6  Given  the  preconditions  of  Lemma  3.f,  polarr(r^)  C  ,  where  r  <  m.m.^^/s/2,  S/5). 

Proof:  Recall  that  polarr(r^)  =  finer5  polarr(u)  and  =  n„er6  hv,  where  hv  was  defined  earlier.  For 

each  !)  £  rj,  polarr(u)  and  hv  are  halfspaces  with  parallel  supporting  hyperplanes  (both  on  the  same 
side  of  the  origin),  which  may  differ  only  with  respect  to  their  distances  from  the  origin. 

We  consider  the  two  cases  arising  in  the  definition  of  hv.  If  ||i>||  >  <5/3,  then  hv  is  at  distance  5/6  from 
the  origin  and  polarr(u)  is  at  distance  r2/||u||  <  (<5/5)2/ (<5/3)  <  5/ 6,  which  is  at  least  as  close.  On  the 
other  hand,  if  ||v||  <  5/3,  then  hv  is  at  distance  £/(2||i>||)  and  polarr.(u)  is  at  distance  r2/||u||  <  £/(2||u||), 
which  is  also  at  least  as  close.  Since  each  bounding  halfspace  of  polarr(r^)  is  as  close  to  the  origin  as  its 
parallel  bounding  halfspace  of  it  follows  that  polar,,  (r*)  C  rgd.  □ 

Given  the  above  results,  we  can  now  prove  the  first  part  of  Lemma  3.4  in  a  manner  analogous  to  the 
proof  of  Lemma  3.1. 

Proof:  (of  Lemma  3.4) 
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Recall  that  Let  0  <  £  <  1  and  5  <  3y/e.  Since  r5  is  a  subset  of  the  orthogonal  projection  of  onto 
the  plane  h(q'),  area(r^)  <  area(Z?95).  Let  r  =  <5/5.  Since  <5  <  3 y/e,  r  <  3^/e/5  <  \Jej2.  Thus,  r  satisfies 
the  conditions  of  Lemma  3.6.  By  applying  both  Lemmas  3.5  and  3.6,  we  have 

area(Vor (D^)dH(q))  >  area(r^)  >  area(polarr(r^)). 

As  observed  earlier,  polarr(T^)  is  a  scaled  copy  of  polar(T5)  by  a  factor  of  r2  =  (5/ 5)2.  Since  this  is  a 
(d—  l)-dimensional  set,  we  have  area(polar,,(r5))  =  (<5/5)2^-1^  •  area(polar(T5)).  As  before,  applying 
the  inequality  from  (16)  on  the  Mahler  volume,  there  is  a  constant  c!  depending  only  on  d  such  that 
area(r^)  •  area(polar(r^))  >  c' .  Combining  these  observations,  we  have 


area(Ugd)  •  area(Vor(Dg'5)  ft  H(q)) 


>  area(r^)  ^ 


2(d— 1) 


area(polar(T5)) 


> 


52(d-l) 


•  J2(d_1). 


Setting  c  =  d / (52^d  1-))  completes  the  proof  of  Lemma  3.4(i).  The  second  part  is  established  by  the  same 
argument  used  in  the  proof  of  Lemma  3.1(h).  □ 


4  Polytope  Approximation 


In  this  section  we  consider  how  the  results  of  the  previous  section5  can  be  applied  to  polytope  approx¬ 
imation.  Our  main  result  is  presented  in  the  following  theorem,  which  establishes  a  uniform  bound 
on  the  complexity  of  convex  approximation  based  on  area,  rather  than  diameter.  As  mentioned  in  the 
introduction,  this  bound  is  never  significantly  worse  that  Dudley’s,  but  may  be  significantly  better  for 
skinny  convex  bodies.  The  proof  appears  later  in  this  section. 


Theorem  4.1  Let  K  C  R d  be  a  convex  body  such  that  the  width  of  K  in  any  direction  is  at  least  e. 
There  exists  an  e- approximating  poly  tope  P  with  number  of  facets 


O(rlogr),  for  r 


i/area(A') 
£•0-1)72  ' 


(This  can  be  adapted  to  show  that  there  exists  an  approximating  polytope  with  O(rlogr)  vertices 
instead  of  facets.)  Note  that  given  that  an  approximating  polytope  P  with  O(rlogr)  facets  exists,  then 
a  polytope  with  at  most  0(rlog2r)  facets  can  be  constructed  using  (9).  In  order  to  prove  the  above 
theorem,  we  start  by  stating  a  lemma  that  connects  dual  caps  and  polytope  approximation.  The  proof 
follows  from  the  works  of  Dudley  (12)  and  Bronshteyn  and  Ivanov  (5).  An  analogous  lemma  relates  caps 
and  poly  tope  approximation. 


Lemma  4.2  Let  N  be  any  set  of  points  on  dK  that  stabs  the  set  of  all  e-dual  caps  of  K.  Then  the 
polytope  defined  by  tangent  planes  constructed  at  the  points  of  N  is  an  e- approximation  of  K. 


In  order  to  use  £-nets  to  stab  all  dual  caps,  we  need  to  show  that  the  range  spaces  of  dual  caps  have 
constant  VC-dimension.  Throughout,  S  denotes  the  sphere  of  radius  3  centered  at  the  origin. 

Lemma  4.3  Let  K  C  Rd  be  a  convex  body  contained  within  the  unit  ball  centered  at  the  origin.  All  the 
following  range  spaces  (X,R)  have  VC-dimension  at  most  d+1. 

5  The  results  of  Section  3  rely  on  the  assumption  that  K  is  e-smooth.  Given  a  general  (non-smooth)  convex  body,  we 
can  apply  a  simple  reduction  to  the  smooth  case,  as  follows.  First,  we  dilate  K  by  performing  a  Minkowski  sum  with  a 
Euclidean  ball  of  radius  e/2.  Let  K'  be  the  result.  Observe  that  the  Hausdorff  distance  between  K'  and  I\  is  e/2  and  K'  is 
(e/2)-smooth.  Under  our  assumption  that  the  width  of  K  in  any  direction  is  at  least  e,  this  dilation  does  not  increase  the 
surface  area  of  K  by  more  than  a  constant  factor.  Next,  we  apply  the  results  of  this  section  and  subsequent  approximations, 
but  with  the  parameter  e  set  to  e/2.  The  combination  of  the  two  error  bounds  (e/2  from  dilation  and  e/2  from  subsequent 
approximations)  results  in  a  total  error  of  e.  The  modification  of  the  e  term  only  affects  the  constant  factors. 
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(1)  X  =  dK  and  R  is  the  set  of  e-dual  caps. 

(2)  X  =  S  and  R  is  the  set  of  Voronoi  cells  of  the  e-dual  caps. 

Proof:  First,  we  consider  an  additional  range  space  define  by  X  =  dK  and  R  is  the  set  of  e-caps.  An  e- 
cap  only  contains  the  points  that  the  corresponding  halfspace  contains,  and  therefore  has  VC-dimension 
upper  bounded  by  the  VC-dimension  of  d-dimensional  halfspaces,  which  is  d  +  1.  The  range  space  (1) 
is  then  equivalent  to  d-dimensional  halfspaces  by  point-hyperplane  duality.  To  prove  the  lemma  for  (2) , 
note  that  the  Voronoi  mapping  defines  a  bijection  between  dK  and  S,  and  therefore  this  range  space  is 
combinatorially  identical  to  the  previous  one.  □ 

Proof:  (of  Theorem  4.1) 

Without  loss  of  generality,  we  may  assume  that  K  has  been  uniformly  scaled  to  lie  within  a  unit 
ball  centered  at  the  origin  and  that  e  has  been  scaled  accordingly  (because  the  ratio  area(A')/e^~b  js 
unaffected  by  uniform  scaling).  By  Lemma  4.2,  if  we  can  stab  all  e-dual  caps  of  K  with  O(rlogr)  points, 
then  we  have  an  e-approximating  polytope  with  the  same  number  of  facets.  We  distinguish  between  two 
types  of  e-dual  caps.  A  dual  cap  D  is  of  type  1,  if  area(I?)  >  y/area (K)  ■  e^d_1^2,  and  otherwise  it  is  of 
type  2. 

We  use  different  strategies  for  stabbing  the  two  types  of  dual  caps.  For  dual  caps  of  type  1,  we  use 
standard  machinery  to  build  a  (l/r)-net  Ni  of  size  0(r  log  r)  in  the  range  space  (1)  of  Lemma  4.3.  The 
net  N\  stabs  all  dual  caps  of  type  1  because  of  their  sufficiently  large  area. 

From  Lemma  3.1,  for  a  type-2  dual  cap  D  we  have  area(Vor(Z?)  fl  S)  =  n(e^_1^2/y/ area(AT)).  In 
order  to  stab  type-2  dual  caps,  we  build  an  0(l/r)-net  V2  for  the  range  space  (2)  of  Lemma  4.3.  The 
net  N2  stabs  the  Voronoi  cells  of  all  type-2  dual  caps.  We  then  select  the  nearest  neighbor  on  dK  for 
each  point  of  N2,  obtaining  a  set  N2.  The  set  N2  stabs  all  dual  caps  of  type  2,  which  concludes  the 
proof.  □ 


5  Polytope  Membership  Queries 

Consider  a  convex  polytope  K  C  [— l,l]d.  In  polytope  membership  queries,  we  are  given  a  query  point 
q,  and  are  to  determine  whether  q  lies  inside  K  subject  to  an  error  bound  e.  Points  within  distance  e  of 
the  boundary  of  K  may  be  reported  as  either  inside  or  outside  of  K .  A  simple  and  natural  algorithm, 
called  SplitReduce,  was  presented  in  (2)  for  preprocessing  K  in  order  to  answer  such  queries  efficiently. 
Although  an  initial  analysis  was  given  in  (2),  the  complexity  of  the  resulting  data  structure  remains  a 
challenging  problem.  In  this  section,  we  show  that,  without  increasing  storage,  it  is  possible  to  improve 
the  query  time  bounds  to  essentially  the  square  root  of  those  given  in  (2) . 

Before  describing  the  algorithm,  we  need  a  definition.  We  say  that  P  is  an  ^-approximation  of  K 
within  a  quadtree  box  Q  if,  for  any  point  x  £  Q,  the  following  holds:  (a)  if  x  £  K,  then  x  £  P,  and  (b) 
if  x  is  at  distance  at  least  £  from  K,  then  x  qL  P.  Note  that  this  definition  implies  that,  for  any  query 
point  x  £  Q,  we  can  correctly  answer  £-approximate  polytope  membership  queries  with  respect  to  K  by 
checking  whether  x  £  P. 

The  SplitReduce  algorithm  is  presented  below.  It  is  given  the  polytope  K  C  [—1,  l]d,  the  approxima¬ 
tion  parameter  e,  and  the  desired  query  time  t.  The  initial  quadtree  box  is  Q  =  [—1,  l]d. 

SplitReduce(Q): 

1.  Let  P  be  an  ^approximation  of  K  within  Q. 

2.  If  the  number  of  facets  |P|  <  t,  then  Q  stores  the  hyperplanes  bounding  P. 

3.  Otherwise,  split  Q  into  2d  quadtree  boxes  and  invoke  SplitReduce  on  each  such  box. 

The  output  of  SplitReduce(Q)  is  a  quadtree  whose  leaf  cells  induce  a  subdivision  of  Q.  Each  leaf  cell 
Q  stores  a  set  of  halfspaces  of  size  at  most  t  whose  intersection  approximates  K  within  Q.  Let  t(Q) 
denote  the  size  of  this  set.  The  storage  of  this  quadtree  is  defined  as  the  total  number  of  halfspaces 
stored  over  all  the  leaf  cells.  A  query  is  answered  by  performing  a  point  location  in  the  quadtree  followed 
by  a  brute  force  inspection  of  the  halfspaces  in  the  resulting  leaf  node.  Because  cells  of  diameter  e  are 
not  subdivided,  the  height  of  the  quadtree,  and  hence  the  time  for  point  location,  is  0(log(l/£)).  For 


12 


the  sake  of  analysis,  let  us  make  the  simplifying  assumption  that  the  number  t(Q)  of  facets  computed  in 
Step  1  is  minimum.  (If  instead,  Clarkson’s  approximation  algorithm  is  used  (9),  the  query  time  increases 
by  a  negligible  factor  of  0(log(l/e)).) 

The  main  result  of  this  section  is  the  following  space-time  tradeoff  for  SplitReduce. 

Theorem  5.1  Given  a  polytope  K  of  diameter  1,  an  approximation  bound  e  >  0  and  a  real  constant  a  > 
2,  it  is  possible  to  answer  e- approximate  polytope  membership  queries  in  time  t  =  l/£(d_1V(2“)  log(l/e)  > 
1  from  a  data  structure  of  space 

O  ^(1/e)  (d~1}  (1_  2L>°S2  4#^  1 ) 

In  order  to  establish  the  tradeoff,  we  first  show  how  to  generalize  Lemma  4.2  to  the  technically  more 
challenging  context  of  approximating  K  within  a  box  Q. 

Lemma  5.2  Let  K  be  a  convex  body  (of  arbitrary  diameter).  Let  Q  be  a  quadtree  box  that  intersects 
dK.  Let  E  be  the  portion  of  dK  that  is  within  distance  yfe  of  Q.  Let  N  be  any  set  of  points  that  stab  the 
set  of  all  y/E -restricted  e-dual  caps  contained  within  E,  and  let  P  be  the  polytope  defined  by  the  tangent 
planes  constructed  at  the  points  of  N.  Then  P  is  an  e- approximation  of  K  within  Q. 

Proof:  Let  x  be  any  point  in  Q.  Since  P  is  defined  by  the  planes  tangent  to  K,  it  follows  that  ii  x  £  K 
then  x  £  P.  Next  suppose  that  x  is  at  distance  at  least  e  from  K.  We  will  show  that  x  (jL  P,  which 
will  prove  the  lemma.  Consider  the  line  segment  joining  x  to  the  point  in  K  n  Q  that  is  closest  to  it. 
By  continuity,  there  must  be  a  point  q  on  this  segment  that  is  at  distance  exactly  e  from  K.  Since  Q  is 
convex,  this  segment  and,  hence,  q  must  be  contained  in  Q.  Consider  the  -^/e-restricted  dual  cap  whose 
apex  is  q.  Note  that  this  dual  cap  is  contained  within  E,  and  so  there  is  a  point  of  N  that  stabs  it. 
Clearly,  the  tangent  plane  at  this  point  separates  q  from  K ,  and  hence  x  from  K.  Thus  x  fL  P,  which 
completes  the  proof.  □ 


As  in  (2),  we  start  by  analyzing  the  low  space  case.  Theorem  5.1  can  then  be  proved  by  induction. 
Our  key  idea  is  to  use  the  machinery  we  developed  in  Section  3.  The  proof  of  the  following  lemma  is 
similar  to  the  proof  of  Theorem  4.1,  but  using  Lemma  3.4  (with  5  set  to  yje)  instead  of  Lemma  3.1. 
Recall  that  S  is  the  Dudley  hypersphere,  which  we  defined  as  the  (d  —  l)-dimensional  sphere  of  radius  3 
centered  at  the  origin. 


Lemma  5.3  Let  K,  Q  and  E  be  as  in  Lemma  5.2,  such  that  K  is  e-smooth  and  Q  is  of  diameter  at  least 
\fe.  Further,  let  E  be  contained  within  a  ball  of  unit  radius  centered  at  the  origin.  Then  it  is  possible  to 
e-approximate  K  within  Q  by  a  poly  tope  with  number  of  facets 


O(rlogr),  for  r 


i/area (E)  •  area(Vor (E)  n  S) 
£(d- 1)/2 


Proof:  By  Lemma  5.2,  if  we  can  stab  the  set  of  all  i/e-restricted  e-dual  caps  contained  within  E 
with  0{r  log  r)  points,  then  we  have  an  e-approximating  polytope  within  Q  with  0(r  log  r)  facets.  To 
prove  this  lemma,  we  distinguish  between  two  types  of  such  dual  caps.  A  dual  cap  D  is  of  type  1,  if 
area(D)  >  y/area(S)/area(Vor(Ll)  n  S)  ■  e^-1)/2,  and  otherwise  it  is  of  type  2. 

By  a  straightforward  generalization  of  Lemma  4.3,  it  follows  that  both  types  of  restricted  dual  caps 
are  of  constant  VC-dimension.  We  use  different  strategies  for  stabbing  the  two  types  of  dual  caps.  For 
dual  caps  of  type  1,  we  use  standard  machinery  to  build  a  (l/r)-net  Ni,  of  size  O(rlogr)  for  the  range 
space  (1)  of  Lemma  4.3,  considering  only  the  portion  E  of  dK.  The  net  stabs  all  dual  caps  of  type 
1  because  of  their  sufficiently  large  area. 

It  follows  from  Lemma  3.4  that  for  a  type-2  dual  cap  D ,  we  have 


area  (Vor(D)  nS)  =  fi 


r(d— 1)/2  area  (Vot(E)  n  S) 
y  area  (E) 


In  order  to  stab  type-2  caps,  we  build  an  0(l/r)-net  N2  of  size  O(rlogr)  for  the  range  space  (2)  of 
Lemma  4.3,  considering  only  Vor (E)  n  S.  The  net  N2  stabs  the  Voronoi  cell  of  all  type-2  caps.  We  then 


13 


select  the  nearest  neighbor  on  OK  for  each  point  of  iV2,  obtaining  a  set  N2.  The  set  N2  stabs  all  caps  of 
type  2,  which  concludes  the  proof.  □ 


The  following  lemma  will  be  useful  in  establishing  the  storage  bounds.  Throughout  the  remainder  of 
the  analysis,  Q  will  denote  a  quadtree  box,  and  T  will  denote  the  quadtree  produced  by  SplitReduce(Q) 
for  some  choice  of  parameters  t  >  1  and  e  <  diam(Q).  Given  T,  let  T'  denote  the  subtree  of  T  induced 
by  cells  of  diameter  at  least  e  diam(Q)/2.  Let  L\  denote  the  leaf  cells  of  T'  that  are  not  subdivided 
further  (i.e.,  they  are  leaf  cells  of  T),  and  let  L2  denote  the  remaining  leaf  cells  of  T' .  Observe  that  if  Q 
is  any  cell  of  L2,  then  the  diameter  of  Q  is  at  most  y/e  diam(Q),  since  otherwise  its  children  would  have 
diameter  more  than  \J e  diam(Q)/2,  and  Q  would  not  be  a  leaf  cell  of  T' . 


Lemma  5.4 

(i)  Eq6Li^(Q)  =  0((diam(Q)/£)(d-1)/2). 

(a)  \l2\  =  o 

Proof:  The  proof  of  (i)  follows  from  the  proof  of  Lemma  2.2  in  (2).  For  completeness,  we  provide 
a  short  sketch.  Dudley’s  approximation  yields  an  e-approximation  to  I\  n  Q  as  the  intersection  of 
0((diam(<5)/e)^d_1^2)  halfspaces.  In  Dudley’s  construction,  each  of  these  halfspaces  is  responsible  for 
approximating  a  region  of  K’s  boundary  of  diameter  0(y/e  diam(Q)).  By  a  simple  packing  argument, 
each  halfspace  is  useful  in  this  manner  for  only  a  constant  number  of  cells  Q  €  L 1,  which  establishes 
assertion  (i). 

To  prove  (ii),  let  Q  be  any  cell  of  L2 .  Recall  that  ^/e  diam(Q)/2  <  diam(Q)  <  y/e  diam(Q).  Let 
e  =  e/(4  diam(Q)).  For  the  remainder  of  this  proof,  we  take  the  center  of  Q  as  the  origin  and  assume 
that  we  have  scaled  space  so  that  Q  has  diameter  1.  Note  that  in  this  scaled  space  <  diam(Q)  <  2 
As  usual,  we  let  S  denote  the  Dudley  hypersphere,  defined  as  the  (d  —  l)-dimensional  sphere  of  radius  3 
centered  at  the  origin.  Let  K  be  the  Minkowski  sum  of  K  and  a  ball  of  radius  e  (thus,  I\  is  e-smooth). 
Let  E  denote  the  portion  of  dK  that  is  within  distance  at  most  of  Q.  Note  that  E  is  contained 
within  the  unit  ball  centered  at  the  origin. 

The  intuition  underlying  our  approach  for  bounding  the  number  of  cells  in  L2  is  to  show  that 
area(Vor(i?)  f~l  S)  must  be  large.  Since  the  area  of  S'  is  a  constant,  we  will  show  that  this  implies 
that  there  cannot  be  too  many  cells  in  L2.  In  order  to  show  that  area(Vor(i?)  (iS)  is  large,  we  will  apply 
Lemma  5.3  and  show  that  otherwise  we  could  have  approximated  K  within  Q  to  the  desired  accuracy 
with  fewer  than  t  facets  and  Q  would  then  not  have  been  subdivided  further.  We  now  make  these  ideas 
precise. 

We  first  prove  that  area(Vor(A)  fl  S)  =  fi((t/  log  t)2  ■  e^d_1^2).  For  the  sake  of  contradiction,  suppose 
that  area(Vor(i?)  fl  S)  <  c  •  (f/logf)2  •  e^d_1)/2  for  a  suitable  constant  c.  Clearly,  the  conditions  of 
Lemma  5.3  are  satisfied  with  K,Q,E,  and  e  in  place  of  K,Q,E,  and  e,  respectively.  Therefore,  K  can 
be  e-approximated  within  Q  by  a  polytope  with  m  facets,  where 


to  <  Cirlogr,  for  r  = 


y area (E)  ■  area  ^Vor(A)  fl  S^j 


eid- 1)/2 


where  c\  is  the  constant  hidden  in  the  O-notation  of  Lemma  5.3.  Note  that  an  e  absolute  error  in  the 
space  after  scaling  corresponds  to  e-diam(Q)  =  e/4  absolute  error  prior  to  scaling.  Thus,  prior  to  scaling, 
to  facets  suffice  to  provide  an  e-approximation  of  K  within  Q.  To  compute  to,  observe  that 

area (E)  <  c2  diam(Q)d"1  <  c2  2d~l  e< d~1] /2, 


where  c2  is  some  constant.  By  our  assumption,  area(Vor(I?)  fl  S')  <  c  •  (t/logt)2  •  e/d  1 Substituting 
these  upper  bounds  on  areas,  we  obtain 

sjc2  2d~1  £(d-l)/2  .  c  (t/fogt)2  gtd-l)/2 
r  -  e(rf-i)/2 

<  V2 d~1c2c- 

log  t 
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Recalling  that  m  <  cirlogr,  it  is  easy  to  see  that  we  can  make  m  <t  by  choosing  c  to  be  sufficiently 
small.  Therefore,  Q  should  not  have  been  subdivided  further,  which  contradicts  the  fact  that  Q  £  L2. 
Thus  area(Vor(U)  n  S)  =  fl((f/logf)2  •  gf^-1)/2)  =  Sl((f/logf)2  •  (£/diam(<5))(d_1^2). 

Since  the  area  of  S'  is  a  constant  and  any  point  of  S  can  be  contained  in  Vor (E)  for  at  most  a  constant 
number  of  cells  Q  £  L2,  it  follows  that  the  cardinality  of  L2  is  0((logt/<)2  •  (diam(<5)/£)(d_1^2),  which 
completes  the  proof.  □ 


We  are  now  ready  to  present  our  improved  analysis  of  the  low  space  case. 

Lemma  5.5  The  output  of  SplitReduce(Q)  for 

/diamWlV-1— ),8lo  /diam(Q)\  1 

is  a  quadtree  with  storage  0((diam(<3)/£)(d_1^2). 

Proof:  Let  T  denote  the  quadtree  produced  by  the  algorithm.  We  will  show  that  ^f(Q )  over  all  leaf 
cells  Q  of  T  is  0((diam(Q)/£)(d_1)/2),  which  will  establish  the  desired  storage  bound.  Recall  T' ,Li,  and 
L2  as  defined  just  before  Lemma  5.4.  By  Lemma  5.4(i),  J2qeLi  KQ)  =  0((diam(Q) /£)^d_1^2). 

Next,  consider  the  set  L2  of  cells.  Let  Q  be  any  cell  of  L2.  Recall  that  yje  diam(Q)/2  <  diam(Q)  < 
\Je  diam(Q).  Since  t  >  (diam(Q)/£)(d_1)/8,  it  follows  that  t  >  (diam(Q)/£)^d_1^4  >  1  for  £  < 
diam(<5)/4  (if  £  >  diam(Q)/4  then  the  lemma  holds  trivially). 

Lemma  2.2  in  (2)  states  that  the  output  of  SplitReduce(Q)  for  t  >  (diam(Q)/£)^d_1^4  >  1  is  a 
quadtree  with  storage  0((diam(Q)/£)(d_1^/2).  Therefore,  the  storage  for  the  subtree  rooted  at  Q  is  at 
most  0((diam(<5)/£)(d_1^2)  =  0((diam(Q)/£)(d-1)/4).  Also,  by  Lemma  5.4(h),  |L2|  =  0((logf/f)2  • 
(diam((5)/£)^d_1^2)  =  0((diam(Q)/£)^d_1^4),  where  we  have  used  the  lower  bound  on  t  given  in  the 
statement  of  the  lemma.  It  follows  that  the  total  storage  for  the  subtrees  rooted  at  cells  of  L2  is  at  most 
0((diam((2)/£)(d-1)/4  •  (diam(Q)/£)(d_1^4)  =  0((diam((5)/£)^d_1^2).  This  completes  the  proof.  □ 

Using  Lemma  5.5  as  the  base  case,  we  now  extend  the  space-time  tradeoff  to  other  query  times. 

Theorem  5.6  Let  a  >2  be  a  constant.  The  output  of  SplitReduce(Q)  for 

t  >  ^diam(Q)^  {d~1)/(-2a)  lQg  ^diam(Q)^  >  1 

is  a  quadtree  with  storage 

Q  ^  ^ diam(Q)  y ' d~ 1)  (4~  2 u°s2  “  Llog2QTJ  _1 )  ^ 

Proof:  Let  k  =  [log2aJ.  Our  proof  proceeds  by  induction  on  a  constant  number  of  steps  k.  The  base 
case  k  =  1  corresponds  to  Lemma  5.5.  Next,  assume  that  the  theorem  holds  for  1, . . . ,  k  —  1,  that  is,  for 
a  <  2k.  We  need  to  prove  that  the  theorem  holds  for  2 k  <  a  <  2fc+1. 

Let  T  denote  the  quadtree  produced  by  the  algorithm  with  2k  <  a  <  2k+1.  Recall  T\  L\,  and  L2  as 
defined  just  before  Lemma  5.4.  By  Lemma  5.4(i),  J2q£Li  KQ)  =  0((diam(Q)/£)0-1)/2)i  and  thus  the 
storage  for  the  leaf  cells  of  L\  satisfies  the  bound  of  the  theorem. 

We  next  consider  the  set  L2  of  cells.  Let  Q  be  any  cell  of  i2.  Recall  that  y/e  diam(Q) /2  <  diam(Q)  < 
1 Je  diam(Q).  It  follows  that 


t  > 


diam(Q) 


(d-l)/(2c/) 

log 


^diam(Q)^ 


>  1 
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where  a'  =  a/2  for  e  <  diam(Q)/16  (if  e  >  diam(<3)/16,  then  the  theorem  holds  trivially).  Therefore, 
we  can  use  the  induction  hypothesis  to  obtain  the  following  storage  for  each  cell  Q  €  L2: 


0  |  ^diam(Q)  ^ 


=  o 


'diam(Q) 


c  ')(j-av-vr 


Also,  by  Lemma  5.4(h),  \L2\  =  0((logf/f)2  •  (diam(Q)/e)(d_1)/2)  =  0((diam(<3)/ff)^d^1H2_a)))  where 
we  have  used  the  value  of  t  given  in  the  statement  of  the  theorem.  It  follows  that  the  total  storage  for 
the  subtrees  rooted  at  cells  of  L2  is  at  most 


O 


'diam(Q) 


This  completes  the  proof.  □ 

Recalling  that  the  time  for  locating  the  leaf  cell  containing  the  query  point  is  0(log(l/£)),  Theorem  5.1 
follows  as  an  immediate  consequence. 
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